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The Belinskii, Khalatnikov and Lifshitz conjecture [1| posits that on approach to 
a space-like singularity in general relativity the dynamics are well approximated by 
'ignoring spatial derivatives in favor of time derivatives.' In [3] we examined this 
O , idea from within a Hamiltonian framework and provided a new formulation of the 

conjecture in terms of variables well suited to loop quantum gravity. We now present 



' the details of the analytical part of that investigation. While our motivation came 

from quantum considerations, thanks to some of its new features, our formulation 
should be useful also for future analytical and numerical investigations within general 
relativity. 



Q . PACS numbers: 04.20.Dw,04.60.Kz,04.60.Pp,98.80.Qc,04.20.Fy 

I. INTRODUCTION 

> 

]^ ■ Originally formulated in 1970, the Belinskii-Khalatnikov-Lifshitz (BKL) conjecture states 

that as one approaches a space-like singularity, 'terms containing time derivatives in Ein- 
stein's equations dominate over those containing spatial derivatives' [l|. This implies that 
^ \ Einstein's partial differential equations are well approximated by ordinary differential equa- 
O ■ tions (ODEs), whence the dynamics of general relativity effectively become local and oscil- 
latory. The time evolution of fields at each spatial point is well approximated by that in 
homogeneous cosmologies, classified by Bianchi [sf. The simplest of these are the Bianchi 
I metrics which have no spatial curvature and the Bianchi II metrics which have 'minimal' 
spatial curvature. According to the BKL conjecture, the dynamics of each spatial point 
• follow the 'Mixmaster' behavior — a sequence of Bianchi I solutions bridged by Bianchi II 
transitions. Finally, with the significant exception of a scalar field, matter contributions 
become negligible — to quote Wheeler, "matter doesn't matter" . 

In the beginning, the conjecture seemed to be coordinate dependent and rather implau- 
sible. However, subsequent analysis by a large number of authors has shown that it can be 
made precise and by now there is an impressive body of numerical and analytical evidence 
in its support [i']. It is fair to say that we are still quite far from a proof of the conjecture in 
the full theory. But there has been outstanding progress in simpler models. In particular, 
Berger, Garfinkle, Moncrief, Isenberg, Weaver and others showed that, in a class of models, 
as the singularity is approached the solutions to the full Einstein field equations approach 
the 'Velocity Term Dominated' (VTD) ones obtained by neglecting spatial derivatives 048[. 
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Andersson and Rendall [9| showed that for gravity coupled to a massless scalar field or a 
stiff fiuid, for every solution to the VTD equations there exists a solution to the full field 
equations that converges to the VTD solution as the singularity is approached, even in the 
absence of symmetries. These results were generalized to also include p-form gauge fields 



m 



10] . In these VTD models the dynamics are simpler, allowing a precise statement of 
the conjecture that could be proven. In the general case, the strongest evidence to date 
comes from numerical evolutions. Berger and Moncrief began a program to analyze generic 



cosmological singularities While the initial work focused on symmetry reduced cases 



12[, more recently Garfinkle [13| has performed numerical evolution of space-times with no 
symmetries in which, again, the Mixmaster behavior is apparent. Finally, additional support 
for the conjecture has come from a numerical study of the behavior of test fields near the 



singularity of a Schwarzschild black hole [14 



With growing evidence for the BKL conjecture, it is natural to consider its implications 
to quantum gravity. The conjecture predicts a dramatic simplification of general relativity 
near space-like singularities, which are precisely the places where quantum gravity effects 
are expected to dominate. A promising approach to analyze this issue is provided by loop 
quantum cosmology (LQC) [l5| where there are now several indications that the quantum 
gravity effects become important only when curvature or matter density are about a percent 
of the Planck scale. Therefore it is quite possible that, generically, spatial derivatives become 
negligible compared to the time derivatives already when the universe is sufficiently classical. 
In this case a quantization of the effective theory with ODEs, that descends from techniques 
applicable in the full theory, could provide a reliable qualitative picture of quantum gravity 
effects near generic space-like singularities. If, on the other hand, the BKL behavior sets in 
only in the Planck regime, this strategy would not be viable. But since there is no reason to 
trust Einstein's equations in this regime, then the conjecture would also not have a physically 
interesting domain of validity. 

LQC is the result of application of the principles of loop quantum gravity (LQG) (l6- isl 
to symmetry reduced cosmological models. Initial study of the k=0 Friedmann Lemaitre 
Robertson Walker (FLRW) models revealed that the quantum geometry effects underly- 
ing LQG provide a natural mechanism for the resolution of the big bang singularity |l9| . 
Subsequent more complete analysis led to a detailed understanding of the physics in the 
Planck regime and also showed that although these effects are very strong there — capable 
of replacing the big bang with a quantum bounce — they die extremely rapidly so as to re- 
cover general relativity as soon as the curvature falls below Planck scale |20|. These results 
were then extended to include spatial curvature in [2l[ and a cosmological constant in j22j. 
More recent investigations reveal that if matter satisfies a non-dissipative equation of state 
P = P{p), LQC resolves all strong curvature singularities of the FLRW models, including. 



e.g., those of the 'big-rip' or 'sudden death' type [23(1 . Also, it is now known in LQC that 



the Bianchi I and II and IX singularities are resolved |24|-|26 



In view of the BKL conjecture, these results, together with further support from the 
'hybrid' quantization of Gowdy models 27| , suggest that there may well be a general theorem 
to the effect that all space-like singularities of the classical theory are naturally resolved in 
LQG. However, it is difficult to test this idea using the current formulations of the BKL 
conjecture since these approaches are motivated by the theory of partial differential equations 
rather than by Hamiltonian or quantum considerations (See e.g., j28, 2^|. In particular, most 



approaches perform a rescaling of their dynamical variables by dividing by the trace of the 
extrinsic curvature. It is difficult to promote the resulting variables to operators on the LQC 
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Hilbert space. In the analysis presented here, we reformulate the BKL conjecture in a way 
better suited to LQC and explore the resulting system both analytically and numerically. 

In LQG one begins with a first order formalism where the basic canonical variables are a 
density-weighted triad and a spin-connection [iMi]. In section mi we will begin by recalling 
this Hamiltonian formulation of general relativity. In section UTTl we rewrite this theory using 
a set of variables that are motivated by the BKL conjecture. Rather surprisingly, the core 
of this theory can be formulated using (density weighted) fields with only internal indices; 
space-time tensors never feature] To understand the implications of the BKL conjecture 
to LQG, we need to express the conjecture using this Hamiltonian framework. This task 
is carried out in section IIVI We provide a weak and a strong version of the conjecture. 
The key idea is to say that, as one approaches space-like singularities, the exact system 
is well approximated by a truncated system which features only time derivatives. Non- 
triviality of the formulation lies in the choice of variables and specification of how limits 
are taken. Our procedure satisfies a number of stringent requirements. In particular, one 
can either first truncate the Hamiltonian and then obtain the equations of motion or first 
obtain the full equations of motion and then truncate them; the two procedures commute. 
In section |V] we study the truncated Hamiltonian system and explore its dynamics in some 
detail. We show that it exhibits all the known features such as the 'u-map' and spikes. 
Thus, the Hamiltonian framework we were led to by LQG considerations successfully captures 
the Mixmaster dynamics faithfully. Therefore, in addition to providing a viable point of 
departure to analyzing the fate of generic space-like singularities in LQG, it should also be 
useful in analytical and numerical investigations of the BKL conjecture in classical general 
relativity itself. In section IVT] we summarize the main results and comment on their relation 
to those of other works. 

The two appendices contain more technical material. Appendix |A] introduces densities 
in a coordinate-free manner. This notion is important because the basic variables in our 
formulation of the BKL conjectures are scalar densities of weight 1. In the main text, for 
simplicity we have set the shift and the Lagrange multiplier of the Gauss constraint equal 
to zero. Appendix |B] contains the full equations without these restrictions. 



II. PRELIMINARIES 

We will consider space-times of the form = M x where is a compact, oriented 
3-dimensional manifold (without boundary).^ We will formulate general relativity in terms 
of first order variables, the point of departure of LQG Jsof. These consist of pairs of fields 
consisting of a (density weighted) orthonormal triad, and its conjugate momentum /T* 
which on solutions will correspond to extrinsic curvature. The fundamental poisson bracket 
is given by 

{Et{x),Kl{y)} = 5l5t5\x-y) (2.1) 

Herein, early letters, a, 6, c denote spatial indices while denote internal indices which 

take values in so(3) — the Lie algebra of S'0(3). Tildes are used to capture density weights 
of quantities; a tilde above indicates that the quantity transforms as a (tensor) density of 



The restriction on topology is made primarily to avoid having to specify boundary conditions and having 
to keep track of surface terms. There is no conceptual obstruction to removing this restriction (following, 
for example, the Hamiltonian framework underlying LQC). 
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weight 1 and a tilde below will denote a (tensor) density of weight —1. The internal indices 
can be freely raised and lowered using a fixed kinematical metric on so(3). The phase 
space spanned by smooth pairs {E^, KD will be denoted by V. 

These phase space variables are related to their Arnowitt, Deser and Misner (ADM) [sij 
counterparts by 

E^E^.q'^ = Iq^'' (2.2) 
K:e1 = ^K,' (2.3) 

where qab is the metric on the leaf ^M, q its determinant, and Kab the extrinsic curvature of 
^M. In terms of these variables we perform a 3+1 decomposition of space-time to obtain as 
Hamiltonian a sum of constraints with Lagrange multipliers [l^, 32 1: 

H[E,K]= [ -lN~S-lN^Va + K&. (2.4) 

The Lagrange multipliers N, N"", the lapse and shift, are related to the choice of slicing and 
time in the standard fashion, and Aj is related to rotations in the internal space. Phase 

space functions S, Va, and & are the scalar, vector, and Gauss constraints (with density 



weights 2, 1, 1 respectively), given by j30|, l32 



SiE,K) ^-~qn-2^%KlKl (2.5) 
Va{E,K) ^Wy,{Ki^&l) (2.6) 
&{E,K) =e/''E]Ki (2.7) 

Where TZ is the scalar curvature of the metric qab- The overall sign and numerical factors 
in the constraints are chosen so they reduce to the standard ADM constraints upon solving 
the Gauss constraint. TZ can be written in terms of the triad and its inverse or in terms of 
the triad and the connection compatible with the triad, which is defined by 

DaE': + e.,,TiE'' = 0, or Ti = -^Eb, DaE'^ e'^' . (2.8) 

(Note that Da acts only on tensor indices; it treats the internal indices as scalars.) Although 
T\ is determined entirely by Ef for now it is convenient to use all three fields K^^ and 
E'l in our classical analysis: In our formulation of the BKL conjecture and K^^ will be 
the relevant degrees of freedom near the singularity, so it is natural to express the theory in 
terms of them. 

The equations of motion are obtained by taking Poisson brackets with the Hamiltonian 
on the phase space V: 

Et = {E^,H[E,K]} (2.9) 
iir: = {KIH[E,K]}. (2.10) 



V is the phase space underlying LQG. The basic variables (A^, E^) used there are obtained 
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by a simple canonical transformation on V 30|: 

{Et,K)^{A^^,r'K) with A\ = V\ + ^Kl. 



(2.11) 



7 being the Barbero-Immirzi parameter of LQC. (In classical general relativity, space-time 
equations of motion are independent of the value of this real parameter.) For sirnplicity 
of presentation we will introduce our formulation of the BKL conjecture using (i?f , KD 
although it will be clear that our framework can be readily recast in terms of {E^, A^). 



III. VARIABLES MOTIVATED BY THE BKL CONJECTURE 

In order to formulate the BKL conjecture in this system, one needs to specify two things: 
What kind of derivatives are to dominate as one approaches the singularity and what kind 
are to become negligible? And, what are the quantities whose derivatives are to be treated as 
negligible? In this section we first motivate and introduce a set of variables and a derivative 
operator and then use them to formulate the conjecture. The main idea is as follows. The 
accumulated evidence to date suggests that the spatial metric qab becomes degenerate at 
the space-like singularity whence its determinant q vanishes there. (In particular, this is 
borne out in the numerical simulations of solutions with two commuting Killing fields — the 
so-called G2 space-times which include Gowdy models j33|.) We will focus on the class 
of singularities where this occurs. In this case one would expect that if we rescaled fields 
which are ordinarily divergent at the singularity with appropriate powers of g, the rescaled 
quantities would have well defined limits. 

Now, the density weighted triad is obtained by rescaling of the orthonormal triad e^, 
which is divergent at the singularity, by ^/q. In examples, the factor of ^/q not only gives E^ 
a well defined limit, but the limit in fact vanishes. Therefore, contraction by Ef can serve to 
tame fields which would otherwise have been divergent at the singularity. This consideration 
leads us to construct scalar densities by contracting E^ with K"*, and V\. As noted above, 
since contraction with E'l will suppresses the divergence of and F^, the combination is 
expected to remain finite at the singularity. Let us then set 

P,^ ■.= EtKi-ElK%^ (3.1) 
C,' ■.= EtTi-ElT%K (3.2) 

These two fields, Pi^ and Ci^ will turn out to be the relevant variables near the singularity 
in our BKL framework. In particular, we will show below that the constraints of general 
relativity can be expressed in terms of polynomials of these basic variables and their deriva- 
tives. Therefore if the basic variables and their derivatives remain finite at the singularity, 
the constraints will also continue to hold there. Since the Hamiltonian of the theory is a 
linear combination of these constraints, dynamics of the basic variables will meaningfully 
extend to the singularity. 

Beyond the possibility of being bounded at the singularity, an important feature of these 
variables is that they have only internal indices which can be freely raised and lowered 
using the fixed, kinematic, internal metric g*-'; the dynamical metric q"'^ which diverges at 
singularities is not needed. Under diffeomorphisms Pj-' and Ci^ transform as density weighted 
scalars on ^M. Because of this feature, statements about their asymptotic properties can 
be formulated much more easily than would be possible if they were tensor fields. (For a 
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coordinate free introduction to densities, see Appendix 

To illustrate why these variables are likely to be well defined at the singularity, let us 
consider the Bianchi I model. Because of spatial fiatness, we can work in an internal gauge 
in which Ci' = everywhere. What about and Pi^7 In terms of the commonly used 
proper time r, the metric is given by ds"^ = —dr'^ + r^^'dxf and the singularity occurs at 
r = 0. Since = 1, we have g = in the Bianchi I chart. In addition, due to the second 
constraint on the exponents, = 1 whence the density weighted triad vanishes at 

the singularity as t^~p^ and P/ is finite there for each i. 

We further introduce a derivative operator Di defined by the contraction of Da with E^: 

Di:=EtDa. (3.3) 

The expectation is that this contraction will have the effect of suppressing terms containing 
Di as we approach the singularity. Thus, Di will be the spatial derivatives we were seeking 
which, when acting on certain quantities, will be conjectured to be negligible near the 
singularity.^ The variable Pij is related to the momentum P"*(conjugate to the 3-metric) 
in the ADM phase space by qP""^ = E^EjP'^^. Cij encodes information in the Di spatial 
derivatives of the triad Ef: 

= -E\ e^^^ DkEI . (3.4) 

Note that, although the Cij depend on spatial derivatives of the triad and are often sub- 
dominant to Pij, it turns out that they are not always negligible in the approach to the 
singularity. Indeed, this behavior is observed in the truncated system, which is discussed 
in section IIVI It is C^^ rather than the triads themselves that will feature directly in our 
formulation of the conjecture. 

For simplicity of notation, from now on we will drop the tildes. Thus, from now on each 
of £'f , Ci^ , Pi^ , Di carries a density weight 1, while the lapse field N carries a density weight 
— 1. The scalar and the vector constraint functions 5* and Vi, (introduce below) carry density 
weight 2 while the Gauss constraint carries density weight 1. 

By making use of (13. ip and (13.41) . functions of {E1,Kl) and their covariant derivatives 
can be rewritten in terms of {Ef, Ci^ , Pi^) and their Di derivatives. The scalar curvature TZ 
for example can be expressed entirely in terms of Ci^ and its Di derivatives: 

q-jZ = -2e''^D,{C,u) - 4q,,]C[^^l - Q,C^' + Uj^ (3.5) 

Consequently, the constraints can be re-expressed entirely in terms of Ci\Pi^ and their Di 



^ This operator is linear and satisfies the Leibnitz rule. It ignores internal indices (since the action of Da 
is non-trivial only on tensor indices) . However since its action on a function / does not yield the exterior 
derivative df, Di is not a connection. If we were to formally treat as a connection, it would have torsion, 
which is related to C: D^iDj^f — —T^jDkf where T\j — ekiliC^y In what follows, Di often acts on 
scalar densities. This action is given explicitly in Appendix El 
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derivatives (with no direct reference to or even the determinant q of the 3- metric): 

S = 2e^^'=A(C,fe) + 4q,,]C[^^] + Qfi^' - + P,,P^' - ^P' (3.6) 

V, = -2D,P,^ - 2e,MP^ - e^.kCP^' + 2e,,kP''Ci' (3.7) 
Qk ^ e'^^Pj,. (3.8) 

Here we have converted the co- vector index on the vector constraint Va to an internal index 
by contracting it with E°i. Since the is assumed to be non-degenerate away from the 
singularity the constraint Vi defines the same constraint surface as the original vector con- 
straint introduced in fl2.5p . Notice here that the constraint can be easily decomposed into 
those terms that contain the derivative Di and those that don't. 

The equations of motion for i?", Ci\ Pi^ can be written in a similar form. These can 
be obtained using the full Poisson brackets fl2.9p .( pTT0|) or by directly computing Poisson 
brackets of Pi^ and Ci^ with the scalar /Hamiltonian constraint. To streamline the second 
calculation, let us specify the Poisson brackets between i??, Ci^ , and Pi^: 



{Et{x),P,\y)} = {E]ix)5,'' - E^ix)6,') 5{x,y) (3.9) 

{P,^(x), Pk\y)] = {Pk'{x)6/ - P/{x)6k')6{x, y) (3.10) 

{/^■P^^ ^QkiC^'} = J{f,,gMiC''6'' + + e^'^^S'^gkiD^fi,) (3.11) 

{E^ix),C^kiy)} = and {Q^(x), Cfc'(y)} = , (3.12) 

where fij,gij are smooth test scalar fields. The equations of motion obtained by taking 
Poisson brackets with the scalar constraint are then given by 

(jij = _ e^f'^DkiNi^SiP - Pi ')) + A^[2C^;pl^'l-'') + 2C^^^^PC - PC'^\ (3.13) 

pn = e^^^Dk{N{l/25\C - Ci')) - e''^"'D^{NCki)S'^ + 2e^''"'C^'''^D^{N) 

{D'D'J - D^Dk5'^)N + N[-2C^'^'^C^^ + CC'^ - 2C^^^^C[ki]5'^] (3.14) 



and 

4« = -NP^'E] 

where we have set the shift to zero to reduce clutter. (For non-zero shift, see Appendix 
[BI ) Note that the equation of motion for E^ is a simple ODE. Note also that, as was the 
case with constraints, the equations of motion for Ci^ and Pi^ can again be written in terms 
of scalar densities and the derivative Di only. This motivates us to ask for an evolution 
equation for the derivative operator Di. Since Di ignores internal indices, it suffices to 
consider its action just on scalar densities 5'„ of weight n. We have: 

71 

DA = - [D,{NP)]Sn - NP,W,Sn . (3.15) 

Thus we have cast all the constraint as well as evolution equations as a closed system 
involving only Ci^, Pi\ and Di. These equations can then be used as follows. On an 
initial slice, we construct (Cj-^, Pi^ , D^) from a pair {E^, KD of canonical variables. Then we 
can deal exclusively with the triplet {Ci^ , Pi^ , Di). The pair {E'^,Kl) satisfies constraints 
if and only if the triplet satisfies fl3.6p - fl3.8p . Given such a triplet, we can evolve it using 
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f l3.13p .( l3.14p . fl3.15p without having to refer back to the original canonical pair (E^jKl). 
These two sets of equations have some interesting unforeseen features. First, as aheady 
mentioned, the basic triplet {Ci^ , Pi^ , Di) has only internal indices: our basic fields are 
scalars on (with density weight 1). It would be of considerable interest to investigate if 
this fact provides new insights into the dynamics of 3+1 dimensional gravity [IBj. Second, 
these equations do not refer to the triad Ef. Suppose we begin at an initial time where 
Ci^ is derived from an E^. Then these constraint and evolution equations ensure that 
Ci^ is derivable from a triad at all times. Furthermore, we can easily construct that triad 
directly from a solution {Ci^ , P/) to these equations: first solve fl3.13p - fl3.15p and then simply 
integrate the ODE 

= -NPi^E^ (3.16) 

at the end. Third, the structure of the constraint and evolution equations in terms of 
{Ci^ , Pi\ Di) is remarkably simple since only low order polynomials of these variables are 
involved. Finally, thanks to our rescaling by ^/q, our basic triplet Ci' ,Pi\Di (as well as 
E^) are expected to have a well behaved limit at the singularity. A close examination of our 
equations shows that they allow the triad become to become degenerate during evolution. 



So, strictly (as in LQG [30|, |32|) we have a generalization of Einstein's equations. 

To summarize, we have found variables which remain finite at the singularity in examples 
and rewritten Einstein's equations as a closed system of differential equations in terms of 
them. Therefore, this formulation may be useful for proving global existence and uniqueness 
results and rigorous exploration of fields near space-like singularities. Finally, although for 
simplicity we have set shift A^* and the smearing field Aj equal to zero, the features we just 
discussed hold more generally (see Appendix |B|) . 

To conclude, let us examine the action of the vector and the Gauss constraints on our 
basic variables. (The action of the scalar constraint yields the evolution equations which 
we have already discussed.) Since the vector constraint generates a combination of spatial 
diffeomorphisms and internal rotations, it is standard to subtract a multiple of the Gauss 
constraint to define the diffeomorphism constraint: 

V: = V,-2{C,^ -^5,^)G, (3.17) 
We can then smear both constraints to obtain 

G[K\ = [ A'^Gk and V'[N] = [ N%' (3.18) 

where A^* s a scalar with density weight —1 so that A^" := N^E^ is the standard lapse and, 
as before. A* has density weight zero. The action of G[A] on the basic variables is given as 
usual via Poisson brackets: 

{P^j,G[A]} = ekijA'Pi' + ekuA'P'', (3.19) 
{C,„ G[A]} = eu.A'Ct + eu^A'C\ + D,A, - Dj^AH,, (3.20) 
{A5'„,G[A]} = e,fe,A^D^S'„ (3.21) 

In the last equation Sn is any scalar density of weight n. As expected the Gauss constraint 
generates infinitesimal 5*0(3) transformations with DiSn and Pj-' transforming as tensors 
and Ci^ transforming as (the contraction of a triad with) an 5*0(3) connection. 
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Similarly, the action of the diffeomorphism constraint is given by the Poisson brackets: 

{P,\ V'[N]] = -2{N'DkP^ + P^W^N^ + P^^eui^N'^C'^) = -2C^P,^ (3.22) 

{C/, V'[N]} = -2{N^DkC,^ + C^i^D^^N^ + Q-''e,,„iV'^C'") = -2£^C/ (3.23) 

{A^n, V'[N]} = -2{NWj{DiSn) + nDj{N^)DA - nt.kiN' C""' D.Sn) = -2C^D,Sn 

(3.24) 

— * 

where N = N"" = E°;Ni. We see that the constraint generates diffeomorphisms as expected 
with Pi^,Ci^, and DiSn transforming as scalar densities. Again, note that the infinitesimal 
changes generated by each constraint involve only the basic variables Ci\ Pi^ , and Di. Thus 
there is still a closed system in terms of this set of variables. 



IV. THE CONJECTURE 

In order to express the BKL conjecture we must make more precise the arena in which 
it is to be applied. The ingredients we need are a space-time with a space-like singularity, 
a notion of 'spatial' and 'temporal' derivatives, and specification of the system to which 
the conjecture is to be applied. We make use of the framework introduced in the previous 
section to provide this arena. 

Let us begin with a 4-manifold, *%f admitting a smooth foliation Mt parameterized by a 
time function, t. We restrict ourselves to a slicing of "^M in which the space-like singularity 
lies on the limiting leaf. This ensures that we can reasonably discuss an approach to the 
singularity as approaching the limiting leaf. The time function t labeling our spatial slices is 
intertwined with the choice of lapse and shift. We will assume that the lapse and the shift 
N\ each with density weight —1, admit a smooth limit as one approaches the singularity. 
Since the spatial metric Qabit) becomes degenerate at the singularity, the commonly used 
lapse function N :— y/qN (with density weight zero) goes to zero, thus placing the singularity 
at t = oo. (These assumptions are minimal and further constraints on admissible foliations 
may well be needed in a more complete framework.) 

Our basic variables will be {Ci^ , Pi-'), the lapse A^, and the shift A^*. By time derivatives, 
we will mean their Lie derivatives along the vector field := TVn" -|- A"" where n" is the 
unit normal to the foliation Mf. By spatial derivatives we will mean their Di derivatives. 
Since Di := E^Da, the notion does not depend on coordinates. Rather, it is tied directly 
to the physical triads and the covariant derivatives compatible with them. Then, the idea 
behind the conjecture is that, as one approaches the singularity, the spatial derivatives 
DiCj'', DiPj^, DiN, DiN^ of the basic fields should become negligible compared to the basic 
fields themselves because of the ^/q multiplier in the definition of which descends to Di. 

We now show that an immediate consequence of this assumption is that the antisymmetric 
part of Cij is negligible compared to the other basic fields. Let us define a* := e^^^Cjk- Then 
by conjecture Di{a^N) is negligible.^ Since the spatial manifold is assumed to be compact. 



^ By their definitions, the internal metric and the alternating tensor Cijk are kinematic, fixed once and 
for all, and are annihilated by all derivative operators Da and 
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integrating this negligible quantity and then integrating by parts, one obtains: 

/ D,{Na')= [ Na'DaE^=l Na'a,, (4.1) 

where we have used the definition of Cij which implies DaE^ = eijkC^^. Since the internal 
metric and the lapse are positive, we conclude that and hence C[jj] are necessarily negligible 
under our assumptions. This fact will be useful throughout our analysis. 

Next, note that we have expressed general relativity in the form of a constrained theory 
in terms of our basic variables, Ci\Pi^ , and Di. Our constraints are composed of quadratic 
terms in our basic variables and terms of the form DiCj^ and DiPj^. We can therefore split 
each constraint into two parts — terms which contain no derivatives and those which do. 
Similarly the equations of motion can be split into terms that contain derivatives and those 
that do not. With this background, we can state two versions of our conjecture. 

Weak Conjecture : As the singularity is approached the terms containing derivatives 
in the constraints and equations of motion are negligible in comparison to the polynomial 
terms. Thus, as the singularity is approached the constraints and equations of motion 
approach those found by setting derivative terms to zero. 

We define the truncated theory to be the system defined by setting Dj-derivative terms 
to zero, 

AC^i' = AP/ = AiV = AiV, = = , (4.2) 

in the equations of motion f l3.13p -( 1XT^ and constraints fl3.6p - fl3.8p . Thus, the weak conjec- 
ture says that the equations of motion can be well approximated by those of the truncated 
theory in the vicinity of the singularity. Note that this does not imply that the solutions of 
the full equations of motion will approach the solutions to the truncated equations as the sin- 
gularity is approached. This additional condition is captured in the strong version as follows. 

Strong Conjecture: As the singularity is approached the constraints and the equations 
of motion approach those of the truncated theory and in addition the solutions to the full 
equations are well approximated by solutions to the truncated equations. 

With the strong conjecture the solution of the full Einstein equations will asymptote to 
solutions of the truncated system defined by (14. 2p . In the following we will analyze this 
truncated system. 

Not only are the truncated constraints purely algebraic, but they involve only quadratic 
combinations our basic variables: 

Q.C^^ - ic^ + p^^p^^ - lp2 (4 3) 

-e,,,CP^' + 2e,jkP''Ci' (4.4) 
e'''Pj^, (4.5) 

The truncated Gauss constraint is in fact exact because (13. 8 p involves no derivative terms, 
while the scalar and the diffeomorphism constraints are genuinely truncated. 

The infinitesimal transformations (I3.19P - (I3.24p generated by the full constraints contain 
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derivative terms that are now assumed to be negligible in comparison to the polynomial 
terms. Ignoring the negligible terms leads us to the following transformations on the basic 
fields: 

{P,„ G(A)}t = 2euuPi)^A' (4.6) 
{Cku G{^)}T = 2ekiuQ)'A' (4.7) 

{P,„ V{N)}t = ^ekiuPi)''N"' (Cj - ^Sj) (4.8) 

{Q„ V{N)}t = 4efc,(,C,)'=iV- (Cj - j5j) (4.9) 

{a„ S{N)}t = -2N (2Cfc(,P^) - PQ,) (4.10) 
{P,„ S{N)}t = -2N i-2akC'j + CQ,) . (4.11) 

The Gauss constraint continues to generate internal rotations, but, whereas in the full 
theory Ci^ transforms as (the contraction of the triad with) a connection, after truncation 
both Ci^ and Pi^ transform as SO (3) tensors. The vector constraint also generates internal 
rotations, since the diffeomorphism constraint generates only negligible terms. 

We arrived at the truncated equations of motion by first obtaining the full equations and 
then applying the truncation to them i.e., by setting spatial derivative terms to zero. But we 
could also have first truncated the constraints to obtain ( 14. 3p - f l4.5p and then computed their 
truncated Poisson brackets with the basic variables. This leads to a consistency check of our 
scheme: do the two procedure yield the same 'truncated equations of motion' in the end? 
The answer is in the affirmative. This fact is illustrated by the following 'commutativity 
diagram': 

Full Constraint Truncation^ Truucatcd Constraint 



Equation of Motion 



Equation of Motion 



Full Equation of Motion "'"'^"'^'^'^^'°'^) Truncated Equation of Motion 



Note that the operation of truncation, the final truncated system, and hence the consis- 
tency requirement mentioned above depends crucially on one's choice of basic variables and 
notions of space and time derivatives. For example, if we had adopted the more 'obvious' 
strategy and used triads E'^ rather than Ci^ as basic variables, we would have been led to set 
Ci^ to zero in the truncation procedure since Ci^ would then be derived quantities, obtained 
by taking the Di derivative of E^. This truncation would have led us just to Bianchi I equa- 
tions. The resulting BKL conjecture would have been manifestly false. Thus, considerable 
care is needed to arrive at variables which satisfy a closed set of equations in a Hamiltonian 
framework, suggest a natural way to make the heuristic idea of ignoring spatial derivatives in 
favor of time derivatives precise, lead to the above commuting diagram, and a version of the 
BKL conjecture that is compatible with the large body of analytical and numerical results 
that has accumulated so far. It is rather striking that the variables {Ci^,Pi^) automatically 
satisfy these rather stringent criteria. 
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V. HAMILTONIAN FORMULATION OF THE TRUNCATED SYSTEM 

In this section we will analyze the truncated system in some detail and show that its 
solutions reproduce the expected BKL behavior. The section is divided into three parts. In 
the first we regard Ci\ Pi^ as fields on the full phase space P, obtain the truncated Poisson 
brackets between them and truncated constraints. In the second we solve and gauge fix the 
vector and the Gauss constraints of the truncated theory. The result is a finite dimensional, 
reduced phase space with a single constraint which is well suited to serve as a starting point 
for quantization inspired by the BKL conjecture. In the third part we discuss several features 
of solutions to this Hamiltonian theory. In particular, we will find that they exhibit Bianchi 
I phases with Bianchi II transitions. 

A. Truncated Poisson brackets 

Since the truncated equations of motion can be formulated entirely in terms of Ci\Pi^, 
let us truncate the Poisson brackets fl3.10p . (13. lip we obtained between them by setting 
the negligible terms on the right side to zero. Since the full Poisson bracket (13. lip involves 
smearing fields fij and Qij we first need to specify which terms involving them are to be 
regarded as negligible. The most natural avenue is to construct to fij and gij only from the 
basic fields (Cj-^, Pi^ , N, Ni, qij, e^^^) (and their Di derivatives). Then the terms containing Di 
derivatives of the smearing fields will also be negligible and hence vanish in the truncation. 
The resulting truncated Poisson brackets between Ci^ and Pi^ are then given by: 

{P,^{x),Ck\y)}T = {Ck'5,' + C^'5,k){x)5{x,y) (5.1) 
{Pi^{x),Pk\y)}T = {Pk'5/-P/6k')ix)6{x,y) (5.2) 
{C^ix),Ck'iy)}T = 0. (5.3) 

These Poisson brackets suffice to determine the equations of motion because the truncated 
Hamiltonian constraint (14. 3 p is algebraic in Ci^ and Pi^ . They are now ODEs, 

C^j = N[2CkiiP''j)-PQj] and Pij = N[-2QkC''j + CQj] , (5.4) 

so the truncated dynamics at any one spatial point decouple from those at other points. 

This system has some notable features. First, we have a closed system expressed entirely 
in terms of Ci^{x) and Pi^{x) at any fixed point x. Furthermore, the equations of motion 
(15. 4 p and constraints (I4.3p -( l43|) are at most quadratic in these variables. In the full theory, 
the triad does not appear explicitly in the equations (I3.13P - (I3.15P but is implicitly present 
through Di. Upon truncation, even this implicit dependence disappears. Second, as in the 
full theory, one can first solve the equations of motion for Ci^{x) and Pi-'{x) and then evolve 
the triad at that point at the end by solving an ODE. Third, the truncated scalar constraint 
(14. 3 p is symmetric under interchange of Ci^ and Pi^ and, by adding a multiple of the Gauss 
constraint, the vector constraint can be made anti-symmetric under this interchange: 

V^T) ■.= e'^''P/CM. (5.5) 

However, this symmetry is broken at the level of equations of motion because the truncated 
Poisson algebra does not have a simple transformation property under this interchange. 
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Because fields at distinct points decouple, to study the truncated system from the view- 
point of differential equations, one can simply restrict oneself to a single spatial point. 
However, this is not directly possible in the Hamiltonian framework because even in the 
truncated theory, the Poisson brackets fl5.ll) . fl5.2p involve 6{x,y). But one can introduce a 
subspace Vhom of the full phase space V tailored to our truncation. Given a point , K^) 
in V consider the pair {Ci^,Pi^) of density weighted fields it determines. The phase space 
point will be said to be homogeneous if there exists an internal gauge and a nowhere van- 
ishing scalar density S^i of weight —1 such that the (density weight zero) scalar fields 
{S^iCi^ , S^iP/) are constants on (and C[ij] = 0). (Fixing a 5*1 is equivalent to fixing a 
3-form on ^M; see Appendix [Al) Clearly, the truncated dynamics leaves this homogeneous 
sub-space Vhom of the phase space invariant. More importantly, Phom is invariant under full 
dynamics: If the Dj-derivatives are initially zero they remain zero under the full equations 
of motion. The Hamiltonian dynamics on "Phom fully captures the truncated dynamics at 
any fixed spatial point on ^M. 

Remark: Since the triads in the full phase space V have been assumed to be non- 
degenerate, they are also non-degenerate in Phom- However, as examples suggest, one would 
expect them to be become degenerate in the limit to the space-like singularity where, how- 
ever, Ci\Pi^ would continue to be well behaved (and some of them may even vanish). It 
is therefore of some interest to extend the homogeneous subspace by adding 'limit points' 
which have this behavior. This construction is not needed in our analysis. However, since it 
may be useful in future investigations, we will conclude this subsection with a brief summary. 
Let us allow the density weighted triads E°; to become degenerate such that the subspaces 
spanned by the non-degenerate directions of vector fields S^iE^ are integrable. (If this con- 
dition is satisfied for one nowhere vanishing scalar density S^i, it is satisfied for all.) Thus, 
in the degenerate case we obtain preferred 2 or 1 dimensional sub-manifolds on ^M. We can 
extend the phase space by including such degenerate E^ if, in addition, the resulting pair 
(Cj-', Pi^) is regular, Cij is symmetric and the pair S-iCi^ , S-iPi^ is homogeneous along the 
preferred lower dimensional sub-manifolds of ^M. Key questions for the BKL conjecture are 
then: i) Does the Hamiltonian flow on V naturally extend to this extension?; and ii) Do 
generic dynamical trajectories flow to it? 

B. Reduced Phase Space 

Since Cij is symmetric but Pij is not, the homogeneous subspace Vhom is not a symplectic 
sub-manifold of the full phase space V. But it turns out that one can obtain a symplectic 
manifold by solving and gauge fixing the truncated vector and the Gauss constraints. It will 
be referred to as the reduced phase space, V^ed- 

The Gauss constraint fl4.5p is equivalent to asking that is symmetric and then the 
vector constraint (14. 4p is equivalent to asking that as matrices, Ci^ and Pi^ should commute. 
To gauge fix the Gauss constraint, we first note the transformation properties (14. 6 p and 
(14. 7p of Pi^ and Ci^ under the action of the Gauss constraint. It is easy to verify that, 
because Pi^ and Ci^ commute, the requirement that they both be diagonal gauge-fixes the 
Gauss constraint completely. It turns out that the diagonality requirement also fixes the 
vector constraint. This may seem surprising at first. But note that the combination V of 
the vector and the Gauss constraint of Eq (15. 5 p again generates internal gauge rotations, 
where, however the generator A* is a 'q-number', i.e., depends on the phase space variables: 
A* = N^{C^j — C6j), where A^-' is the shift used to smear the vector constraint. The fact that 
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the gauge fixing of tlie vector constraint does not impose additional requirements on (Cjj, Pij) 
'cures' tlie mismatcli in tlie degrees of freedom in the homogeneous subspace (arising from 
the fact that while Cij is symmetric, P*-' is not.). 

So far Ci^,P- are fields on ^M, each carrying density weight 1. Since these fields are 
homogeneous, symmetric and diagonal, the reduced phase space is 6 dimensional. It is 
convenient to coordinatize it with just six numbers, Cj, , with J = 1, 2, 3: 

Ci := / Ci^; P^ := / Pi^ etc (5.6) 

where the integrals are well defined because we have completely fixed the internal gauge, in 
that gauge the integrands are all densities of weight 1, and is compact. From now on we 
will focus on the description of Vred in terms of Cj and P^. 

The symplectic structure on Vred is given by the Poisson brackets:"^ 

{P^,pJ} = {Cj,Cj} = and {P',Cj} = 26'jCj. (5.7) 
The scalar or Hamiltonian constraint 

_ CjC^ + ip2 - PjP^ = (5.8) 

now generates the equations of motion via Poisson brackets: 

Pj = NCi{C -2Ci) (5.9) 
Ci = -NCi{P -2Pj) (5.10) 

Here and in what follows we use the summation convention also for the indices J, J and have 
set 

P = Pi + P2 + P3 and C = Ci + C2 + C3. (5.11) 

As a side remark, we note that C/ = is a fixed point of our system for each Cj, whence the 
sign of each Cj along any dynamical trajectory is fixed by the initial conditions. Therefore, 
away from the 'planes' Cj = 0, we can, if we wish, perform a change of variables to Xj = 
In \Ci\/2 and work with the canonically conjugate pair (X^,P/). However, in what follows, 
we will continue to work with (C/, P^). 

Finally, recall that in the BKL conjecture 'the only matter that matters' is a scalar field. 
Let us therefore extend our gravitational reduced phase space to include a massless scalar 
field (f). Denote the conjugate momentum by vr so that {0, tt} = 1. Then on this extended 



^ Note that, thanks to the integrals in the definitions of C/ and , the delta-distributions on the right 
hand side of truncated Poisson brackets (15.11) - (15.21) on V have now disappeared. To write the truncated 
constraints (|4.3I) . (|4.4p in terms of Ci,P^ , one first fixes a nowhere vanishing scalar density 5*1 of weight 1 
(i.e., a 3-form; see Appendix A). One then multiplies these constraints (6*1)^^ to obtain constraints with 
density weight zero. Finally, by noting that Ci = (Ci^5_i) Vo, etc, where Vo is the volume of "^M with 
respect (Si), one obtains the equations of motion for Ci,P^ given below. 
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reduced phase space V^ed the Hamiltonian constraint is given by 

ho^ _ Ci& + _ P,P^ - !^ = (5.12) 

The equations for P and C are still given by (15. 9p and (I5.10p while those of the scalar field 
are simply = tt and tt = 0. 



C. Dynamics 

The Hamiltonian flow in VreA fully captures the gauge invariant properties of the truncated 
dynamics of fields at any one fixed spatial point on ^M. Let us therefore focus on this 
Hamiltonian system. Although the basic constraint and evolution equations on P^ed are 
just ODEs, they have a rich structure; indeed they incorporate the dynamics of all Bianchi 
Type A models. Since the analysis of Bianchi IX is already quite complicated and required 
considerable effort [36|, |37[, we will follow the strategy used in and analyze implications 



of the reduced equations near fixed points. 

There are two sets of fixed points of the dynamics, i.e., points at which Cj = Pj = 0: 

1. Ci = C2, C3 = 0, Pi = P2, P3 = 0, and TT = 

2. Ci = and PjP^ - |p2 + Itt^ = . 

The first set of fixed points corresponds essentially to a dimensional reduction of our theory 



38[ and is therefore highly unstable. To show that our truncation captures the standard 
features associated with the BKL behavior near singularities, it will suffice to focus on the 
second set which, we will now show, in fact corresponds to the Kasner solutions. One can 
show that the solutions to the scalar constraint 2P/P^ — P^ + vr^ = are such that all three 
Pj are positive or all three are negative. Choice of positive signs turns out to be necessary 
and sufficient for the singularity to appear at t = +00 as per our previous conventions. 

Let us return for a moment to the homogeneous phase space Phom and set lapse = S-i, 
the fiducial scalar density for which S-iPij is homogeneous, diagonal, with entries Pi. We 
can then solve the evolution equation (13.161) for the triad -Ef (t) in terms of P/. Finally let 
us set 

2P 

pj = l-^ and r = 6"^*/^ . (5.13) 
Then the space-time metric computed from -Ef (t) is given by 

ds^ = -dr^ + T^P'dxl + T^P^dxj + T^P'dxl (5.14) 
so that the singularity lies at r = (or t = 00). By definition, the constants Pi satisfy 

Pl+P2+P3 = ^ (5.15) 

and the Hamiltonian constraint 



2P/P^ - p2 + = 



(5.16) 
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on P/ translates to the familiar quadratic Kasner constraint 

vl+vl+vl = '^-vl where pi = — . (5.17) 

For each value of < 1, these constraints on the Pi define a 1-parameter family of solutions, 
the intersection of a plane with a 2-sphere. One can check that if p| > 1/2 all the pi are 
positive while if < 1/2 solutions exist only if one of the pi is negative. We will now show 
that this distinction plays the key role for the stability of the solution. 

Let us now move away slightly from a Kasner fixed point (P/, Cj) and consider the 
Hamiltonian trajectory through the new point {Pj, Cj): 

P'j = Pi + 5Pi, C'j = Ci + 5Ci . (5.18) 

Then, the evolution equations for the perturbations are of the form 

{6Pi) = 0{6P^) and {6Ci) = -N6Ci{P - 2Pj) + 0{6C6P) (5.19) 

For definiteness, let us set / = 1. Then P — 2Pi = piP and similarly for / = 2, 3. Now P 
is positive since all three P/ are positive. Therefore, if all Pi are positive (i.e. if p| > 1/2), 

the evolution equation for 6Ci is of the type {SCj) = (negative definite quantity) x 6Ci, 
whence the perturbation will decay, implying stability. In terms of the canonical variables 
describing the scalar field, this occurs when the scalar field is large: Air"^ > P^. This stability 
is in accordance with the Andersson-Rendall results ^ on approach to space-like singularity 
in presence of a massless scalar field in full general relativity. 



CO 
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FIG. 1: Evolution of each of the three Ci in vacuum, starting from a point near the Kasner fixed 
point surface. Initial data is Ci = 1 x 10"^ , C2 = 2 x lO"'^, C3 = 2.2 x lO"'^, Pi = 0.4, P2 = 0.8, 
P3 = 0.0686 (Ci in blue, C2 in green, C3 in red). Since none of the initial C/ vanish, as expected 
from analytical considerations, there is a series of separate Taub transitions between Kasner states. 
Time has been rescaled by a power of 1/4 to allow multiple transitions to be shown on a single 
plot. 

Let us now consider the complementary case where p^ < 1/2. By the above reasoning, 
now (P — 2P/) is negative for some /. For definiteness, let us take Pi to be the largest 
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FIG. 2: Evolution of each of the three Pi in vacuum, starting from a point near the Kasner fixed 
point surface, with the same initial data as in Fig. [1] (Pi in blue, P2 in green, P3 in red). The 
largest eigenvalue, P2, transits first. After this transition. Pi becomes the largest eigenvalue, now 
making Ci unstable. In time all three P/ tend to zero. In terms of parameters pi used in the Kasner 
metric (j5.14p . the initially expanding direction p2 starts contracting at the end of the transition 
and initially contracting pi starts expanding. 



of the P/'s initially so that (P — 2Pi) is negative which implies that Ci will grow and we 
have instability. In this case, we cannot use perturbative analysis for the pair Ci,Pi; it 
is necessary to keep all order terms in Ci and Pi. For simplicity, let us set C2 = C3 = 
initially. Then values of C2, P2,C3, P^ will not change during evolution and equations for 
Ci,Pi simplify, 

P^ = -NCl (5.20) 
Ci = -NCi{P2 + P3 - Pi) = -NCiPpi , (5.21) 

which can be solved exactly to obtain 

Pi(t) = P2 + P3 - 2^P^^ tanh {2^I\P,N{t - Q) 

Ciit) = ±2/^ sech(2/^iV(t - Q) . (5.22) 

These are the Bianchi II solutions written in our variables. Here Ci, the unstable variable, 
rapidly increases and then decays to zero. During that time the Pi transitions between one 
Kasner solution to another. In the asymptotic limits we have: 

Pl(-CX)) = P2 + P3 + 2v/iy^= (v^+ v^)2 

Pi(+oo) = P2 + P3-2v/iV^=(v^-/p;)2. (5.23) 

(In practice the asymptotic limits are achieved quickly, thanks to the hyperbolic functions 
of time.) The result of the transition is that Pi, which was originally was the largest of the 
three P/, has transitioned to a lower value. By a change of variables to the pi used in f l5.14p 
it is apparent that the eigenvalue corresponding to the negative exponent Pi is the one which 
has transitioned, and is positive at the end of the transition. Since the singularity lies at 



18 



r = 0, this means that the initially expanding direction now contracts, and one of the two 
contracting directions now expands. Indeed, fl5.23p is precisely the 'u-map' in pi variables. 

In this analysis we have made the simplification that initially C2 = C3 = 0. If one starts 
from a generic point in the vicinity of the Kasner fixed point set and still with Pi as the 
largest of the three Pj initially, there would again be a transition of the type fl5.22p . But 
as Pi decreases, after a finite time either P2 or P3 will now be the largest eigenvalue and 
making the corresponding Cj unstable. That pair will then evolve according to (15.221) . This 
general scenario was borne out in a large class of simulations of the reduced equations of 
motion. Figs [1] and [2] illustrate this dynamical behavior for generic initial data near the 
Kasner surface. The Taub transitions are easy to see in Fig [2j even though none of the C/ 
are initially zero, the Taub transitions are well described by the analytical expressions (15.221) . 
Fig [3] illustrates the dynamical behavior in cases where the initial data is quite far from the 
Kasner surface. Note that even in this case, the Cj decrease in time so that, although we 
start far away from the Kasner surface, dynamics drives the state to the Kasner surface. 




FIG. 3: Evolution of each of the three Cj in vacuum, starting from a point away from the Kasner 
fixed point surface. Initial data are Ci = 0.7431, C2 = 0.3922, C3 = 0.6555, Pi = 0.1712, 
P2 = 0.0.7060, P3 = —0.3140 (Ci in blue, C2 in green, C3 in red). Even though we start out far 
from the Kasner surface (where all Cj vanish), dynamics drives the state to the Kasner surface. 
Again, time has been rescaled by a power of 1/4 to allow multiple transitions to be shown on a 
single plot. 

We can also draw some lessons for the full theory from this behavior of the truncated 
system. Recall that the dynamical trajectories discussed above can be thought of as rep- 
resenting the evolution of fields at a fixed spatial point. Let us therefore return to and 
consider fields C/(x),P/(x). Now, generically, we will encounter a point xq where the 
all vanish, while being non-zero is the neighborhood of the point. As we noted in section 
IV Bl the sign of C/ is preserved throughout the evolution. This is in particular true during 
the Taub transitions where the magnitude of Cj grows. Therefore, on 'one side' of Xq, a 
Cj will positive and increasing in magnitude, while on the 'other side' it will be negative 
and increasing in magnitude. Therefore its derivative will increase rapidly. Similarly the 
under Taub transitions the values of P/(x) will change from those of one Kasner solution 
to another except at the point xq. Again, this dynamics will generate a large derivative at 
Xq. Thus, analysis of the reduced system suggests that spikes will occur in the full system 
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f l3.13p - f l3.14p . As is well known, these spikes were found in numerical simulations and, more 
recently, also in analytical treatments 43|, |4J| . Whenever spikes appear, the key assumption 
underlying BKL truncation is brought to question because the spatial derivatives are large 
at the spikes. The key issue for the BKL conjecture — and for the application to quantum 
gravity we proposed in section [I] — is whether the time derivatives still dominate generically, 
as they do in examples. 

Let us summarize. Using analytical and numerical methods we showed that there exists 
a well defined subspace "Phom of the full phase space V which exhibits exactly the properties 
expected in the BKL conjecture. Our procedure to arrive at the reduced system is more 



direct than those available in the literature. In [2q], for example, elimination of the off- 
diagonal components and the anti-symmetric parts of Cij involves an additional assumption, 
beyond ignoring spatial derivatives in favor of time derivatives: these quantities are identified 
as part of the 'stable subset' (variables that are expected to decay rapidly as the singularity 
is approached) and then set to zero to obtain the truncated equations. In our treatment, on 
the other hand, the fact that the antisymmetric part of Cij is negligible is directly implied 
by the assumption that the Di derivatives are negligible and the constraints imply that the 
variables Cij and Pij can be simultaneously diagonalized which, furthermore, completely fixes 
the gauge. Thus, our Hamiltonian framework naturally led to a diagonal gauge, enabling 
us to quickly zero-in on the essential variables and eliminating the need to keep track of the 
dynamics of extraneous variables involving frame rotations [28l . 29l |. Finally, the framework 
easily led us to the Mixmaster behavior — a series of Bianchi I phases interspersed by Bianchi 
II transitions. We recovered the 'u-map' for these transitions, and observed the behavior 
expected from the Andersson and Rendall analysis [9| when a scalar field of large enough 
magnitude is introduced. 



VI. DISCUSSION 

We began with the Hamiltonian formulation of general relativity underlying LQG where 
the basic fields are spatial triads with density weight 1, spin connections they deter- 
mine, and extrinsic curvatures K]^. Based on examples that have been studied analytically 
and numerically, it seems reasonable to expect that the determinant q of the spatial metric 
Qab would vanish and the trace K of the extrinsic curvature would diverge at space-like sin- 
gularities. (This expectation is in particular borne out in the numerical simulations of G2 
space-times [ssj.) One can therefore hope to obtain quantities which remain well-defined 
at the singularity either by multiplying the natural geometric fields by suitable powers of q 
or dividing them by suitable powers of K. In the commonly used framework due to Uggla 



et al [13|, |28|, |41|, one chooses to divide by K. One first introduces the so-called Hubble 
normalized triad i^~^e"j by rescaling the orthonormal triad by K~^, and then constructs 
a set of Hubble normalized fields by contracting Fa-' and Kj with K~^e'^i. These fields are 
expected to have regular limit at the space-like singularity. Einstein's equations expressed 
in terms of them naturally suggests a truncation and the truncated system successfully de- 
scribes the expected oscillatory BKL behavior. The resulting form of the BKL conjecture 
is supported by numerical evolutions of full general relativity carried out by Garfinkle (39| . 
However, because there is no underlying Hamiltonian framework, this approach does not 
easily lend itself to non-perturbative quantization. Even if such a framework were to be 
constructed, because of the presence of the factor, it would be difficult to introduce 
quantum operators corresponding to the Hubble rescaled fields. 



20 



Motivated by quantum considerations, we adopted the complementary strategy of multi- 
plying geometrical fields by ^Jq. The LQG Hamiltonian formulation we began with already 
features a density weighted triad with exactly the desired property: E°-i = yge°. Since 
y/q is expected to vanish at the singularity, one can hope to use E"^ in place of the Hubble 
normalized K'^e'} to construct a new set of fields to formulate the BKL conjecture. Indeed, 
(modulo trace terms) our basic variables Ci^ and Pi^ were obtained simply by contacting 
the spatial indices of r^-^ and Kl by E^. Furthermore, because E^ vanishes in the limit, 
the operator Di := EfDa provided a key tool in the formulation of the BKL conjecture: 
asymptotically, DiCj^ and DiPj^ should become 'negligible' relative to Cj^ and Pj^. Now, 
in exact general relativity, time derivatives of Ci^ and Pi^ can be expressed in terms of their 
Di derivatives, purely algebraic (and at most quadratic) combinations of Ci^ and Pj-', the 
lapse and its derivatives (see fl3.6p -f l3T5|) ). Therefore, if in the limit the derivatives 
of the basic fields become negligible compared to the fields themselves, we are naturally 
led to conclude that time derivatives would dominate the spatial derivatives. This chain of 
argument led to our formulation of the BKL conjecture. 

This rather simple idea depends on the fact that the structure of Einstein's equations has 
an interesting and unanticipated feature: as we saw in section UTTl once the triplet Cj-^, Pi' , Di 
is constructed from the triad and the extrinsic curvature i^"* on an initial slice, the 
constraint and evolution equations can be expressed entirely in terms of the triplet. Given a 
solution to these equations, the spatial triad E"; (and hence the metric qab) can be recovered 
at the end simply by solving a total differential equation fl3.16p . This is a surprising and 
potentially deep property of Einstein's equation. It played essential role in our formulation 
of the BKL conjecture and could well capture the primary reason behind the BKL behavior 
observed in examples and numerical simulations. 

Since our framework is developed systematically from a Hamiltonian theory, its BKL 
truncation naturally led to a truncated phase space. The specific truncation used has an 
important property: The truncated constraint and evolution equations on the truncated 
phase space coincide with the truncation of full equations on the full phase space. On 
the truncated phase space we could solve and gauge-fix the Gauss and vector constraints 
to obtain a simple Hamiltonian system (which encompasses all Bianchi type A models). 
Solutions to this system were explored both analytically and numerically. We showed that 
they exhibit the Bianchi I behavior, the Bianchi II transitions and spikes as in the analysis 
of symmetry reduced models j45| and numerical investigations of full general relativity [l3| . 
Therefore, as explained in section [H an appropriate quantization of the truncated system, 
e.g., a la loop quantum cosmology, could go a long way toward understanding the fate of 
generic space-like singularities in quantum gravity. 

In sections UTTl IV Al and IV Bl we restricted ourselves to vacuum equations. The addition 
of a massless scalar field is straightforward and was carried out in the reduced phase space 
framework in section IV Gl If the energy density in the scalar field is small, one again 
has Bianchi II transitions and spikes. However, once the energy density exceeds a critical 
value, these disappear and the asymptotic dynamics at any spatial point is described just 
by the Bianchi I model with a scalar field without transitions. Thus, our truncated system 
faithfully captures the main features generally expected from the analysis of Andersson and 
Rendall [9| in full general relativity coupled to a massless scalar field or stiff fluid. Thus, 
although the initial motivation came from quantum considerations, our formulation of the 
BKL conjecture, and the form of the field equations both in the full and truncated versions, 
should he useful also in the analytical and numerical investigations of singularities in classical 
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general relativity. 

We will conclude with a discussion comparing our approach with that of Uggla, Ellis, 
Wainwright and Elst (UEWE) (jlsl). The Hubble normalized variable used in their formu- 
lation of field equations are give by 

= 3K'^eliKia\j) - K'^elK^Sij (6.1) 

N,, = -3K~'e1J:^a\j) + ^K-'elT% (6.2) 

Ai = -e/^^Sir-^e^r^ (6.3) 

di = 3K-'etda . (6.4) 

These variables are especially useful because they are scale invariant: they are unchanged 
under a constant rescaling of the space-time metric. Because of this property and because of 
the 'regulating' factor in their expressions, it is hoped that in the limit as one approaches 



the space-like singularity, these variables will remain finite j4l|] and their di derivative will 
become negligible. 

We began with quite a different motivation and our focus was on constructing a Hamilto- 
nian framework rather than on differential equations. Since our emphasis was on construct- 
ing phase space variables that can be readily promoted to well-defined quantum operators, 
from the start we avoided the use of factors such as 1/K. As a result, our basic variables 
Ci^ and Pi^ are not scale invariant. Could we have made a different choice which is also well 
suited for quantization and at the same time enjoyed scale invariance? The answer is in the 
negative for the following reason. Under constant conformal rescalings Qab ^^Qab of the 
space-time metric, we have — )■ X^E^, — )■ F^, and K^^ — )■ Kl^. Now, in the analysis of 
approach to singularity, scale invariant quantities are directly useful only if they are space 
scalars and it is not possible to construct scale invariant scalars using just sums of products 
of these fields, i.e., without introducing fields such as for which it is difficult to construct 
quantum operators. Even if one introduces additional non-dynamical fields, such as fiducial 
frames to construct scalars, for natural choices of these frames, scale invariant components 
of fields such as K^, F^ typically diverge at the singularity. Thus, with our motivation, it 
does not seem possible to demand scale invariance of the basic variables that are to feature 
in the BKL conjecture. 

Our viewpoint is that the most important feature of the Hubble normalized variables 
is that although the orthonormal triad typically diverges as one approaches a space- 
like singularity, K diverges even faster, making the combination i^~^e^ go to zero at the 
singularity. Furthermore, it goes to zero at a sufficient rate for its contraction with K^^, F^ 
and da in (16. ip — (16.41) to tame the a priori divergent behavior of these fields. Instead of 
dividing the orthonormal triad e" by K which one expects to diverge at the singularity, 
our strategy was to multiply it by the volume element ^/q which, in examples, goes to zero 
at the singularity. This difference persists also in the treatment of the lapse. The UEWE 
framework assumes that the (scalar) lapse N is such that NK admits a limit N while we 
assume that the density weighted lapse N = {^)~^ N admits a well-defined limit at the 
singularity. Thus, in both cases, the standard scalar lapse N goes to zero so the singularity 
lies at if: = oo. 

The key scale invariant UEWE variables {Nij^T^ij) — which are expected to be well be- 
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haved at the singularity — are related to our {Cij, Kij) via 

= Qp-^ C^ij) and S^^ = -6P"^ P^ij) + 25ij , or, (6.5) 

C{ij) = -^fe) = ~3~ ^^ij ~ (^•^) 

and the two sets of lapse fields are related by: 

^ = K,/^N. (6.7) 

If one focuses only on the structure of differential equations near space-like singularities, 
the two reduced systems would in essence be equivalent if admits a finite, nowhere 

vanishing limit at the singularity. This condition holds for Bianchi I models and also Bianchi 
II which describe the transitions between Bianchi I epochs. In fact in the Bianchi I model, 
^JqK = 1 and our density weighted triad has the same dependence on proper time as the 
Hubble normalized triad. Thus, although the motivations, starting points and procedures 
used in the two frameworks are quite different, surprisingly, in the end the basic variables 
and equations are closely related. 
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Appendix A: Densities 

Since the basic variables that feature in our formulation of the BKL conjecture are scalars 
on of density weight 1, in this Appendix we briefiy recall a coordinate independent 
framework for describe densities. The underlyin g id ea is due to Wheeler and the detailed 
framework was developed by Geroch (see, e.g., |46l|). This framework goes hand in had 
with Penrose's abstract index notation Because the primary application in this 



paper is to our fields Ci^,Pi^ on we will focus on scalar densities on 3- manifolds. But 
generalization to tensor densities on n-manifolds is straightforward. 

Fix an oriented 3-manifold and fix a orientation thereon. Denote by S the space 
of smooth, positively oriented, no-where vanishing, totally skew tensor fields e'**^ on ^M. 
Clearly, given any two elements e"'^'^ and e'"''"^ in S, there exists a (strictly) positive function 
a such that e'"''"^ = ae"'^'^. This fact will be used repeatedly. 

In this paper, a scalar density Sn of weight n is a map from £ to the space of (real valued) 
smooth functions on ^M: e — )■ S'„(e), such that: 



S^{e') = Sn{e) (Al) 

Here n can be any real number but in most applications in general relativity it is an integer. 
(In quantum mechanics, on the other hand, states are (complex-valued) densities of weight 



1/2 on the configuration space j46|.) Since Ci^,Pi^ have density weight 1, let us make a 



short detour to discuss the case n = 1. Fix any 3-form Sabc on ^M. It determines a canonical 



23 



scalar density of weight 1: 

^l(e):=Sa^,ce«'^ (A2) 

Conversely, since 5*1 is a linear mapping from E to smooth functions, it determines a canonical 
3-form Sahc- Thus, our basic variables could also be taken to be 3-forms C'^hc-, P^\hc on 
which take values in second rank tensors in the internal space. The standard ADM 
phase space of general relativity can be similarly coordinatized by positive definite metrics 



qah and tensor fields P^^g which are symmetric in a, h and totally skew in c, (i, e [49|, l50 
Finally note that every metric q^b determines a canonical volume 3-form eahc which has 
positive orientation and satisfies eabc^def q°''^q^^q'^^ = sgn(g) 3!. Therefore it also determines 
a canonical scalar density of weight 1, called the square root of the determinant of qab'- 
^(e) := eabce"'" for all e G £. 

This definition can be extended to density weighted tensor fields in an obvious fashion. 
Note that every carries a natural totally skew tensor density rj"'^'^ of weight 1, called the 
Levi- Ci vita density: 

^-bc^^^^^abc VeG^ (A3) 

Given any metric qab on ^M, the square root of its determinant, ^/q, can also be expressed 
as y/q = T]"^"" eabc- 

Finally, given a derivative operator Da on tensor fields on ^M, we can extend its action 
on densities 5"^ of weight 1 in a natural manner. DaSn is a 1-form with the same density 
weight n, given by 

{DaS^){e) = DaiS^{e))-nXaSn{e) Ve G ^ , (A4) 

where the first term on the right hand side is just the gradient of the function S'„(e) and 
the 1-form Xa is given by Dac'"^'^ = Xac'"^'^. Therefore the action of the derivative operator 
Di introduced in the main text is given by: 

(A5n)(e) = A(5„(e)) - n [E^Xa) S„(e) Ve G ^ . (A5) 

Since the derivative operator Da we considered ignores internal indices, this equation gives 
the action of Di on Ci^ and Pi^ by regarding these basic fields simply as scalar densities with 
weight 1. 

Appendix B: Full Equations of Motion 

In the main text we restricted the equations of motion to the case where the shift is 
zero as is the Lagrange multiplier for the Gauss constraint. In this appendix we give the 
equations of motion in full generality for both full general relativity and in our reduced 
system. The full equations of motion for C and P are as follows. 

C^i = - t^^^Dk{N{]^5\P - Pi ')) + N[2C\P^''^^'> + 2C^^^^PC - PC'^] (Bl) 



+A(A, - N'^Cuj + -N,C) - DkiA" - N'Ci' + -CN'')S,, 
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In the reduced system the derivative terms are set to zero leading to the following equa- 
tions of motion for C and P. 

= TV [2Cfc(,P'=,) - PC,,] + 2eui,C,)\K' - N^Cj + hoN') (B3) 
4 = iV[-2QfcC^- + CQ,] + 2efc;(,P,)'^(A' - N^CJ + \cn') (B4) 
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